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Abstract. We present an extensive experimental study of mode-I, steady, slow crack dynamics in gelatin 
gels. Taking advantage of the sensitivity of the elastic stiffness to gel composition and history we confirm 
and extend the model for fracture of physical hydrogels which we proposed in a previous paper (Nature 
Materials, doi:10.1038/nmatl666 (2006)), which attributes decohesion to the viscoplastic pull-out of the 
network-constituting chains. So, we propose that, in contrast with chemically cross-linked ones, reversible 
gels fracture without chain scission. 

PACS. 62.20.-Mk Mechanical properties of solids - 83. 80. Km Physical gels and microgels - 83. 60. La 
Viscoplasticity, yield stress 

1 Introduction - Physical gels in which cross-linking is due to hydro- 

gen or ionic bonds, much weaker than covalent ones. In 
Hydrogels are a family of materials constituted of a sparse most of them the ne t WO rk is constituted of biopolymers 
random polymer network swollen by a (most often aque- rjj ) e g pro teins (gelatin) or polysaccharides (agar, algi- 
ous) solvent. They can be classified into two subgroups. nates )_ Due to stabilizing steric interactions, these CL may 
- Chemical gels, such as polyacrylamid ones, in which involve many monomeric units (residues), extending over 
the cross-links (hereafter abbreviated as CL) between the lengths of several nanometers. Such is the case for gelatin 
polymer chains are made of single covalent molecular brid- gels. Gelatin results from the denaturation of collagen, 
ges. Their gelation process is irreversible. whose native triple helix structure is locally reconstituted 
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ble segments of single protein chains. Due to the weak et al ^D] have been able to compare the fracture be- 
strength of their CL bonds, physical gels arc thermorc- havior of chemically and physically cross-linked alginate 
versible. For example, gelatin networks "melt" close above gels. They showed that the stiffness/toughness correlation, 
room temperature. This behavior leads to the well studied while agreeing with Tanaka's result for covalent CL, is in- 
slow ageing (strengthening) of their elastic modulus [5], verted for ionic ones. In this latter case "the stiffer the 
and to their noticeable creep under moderate stresses [3]. tougher". 



Biopolymer based physical gels have been attracting 
increasing interest motivated by their wide use in the food 
industry and to promising biomedical developments in 
fields such as drug delivery and tissue engineering All 
these implementations call for the control of their me- 
chanical properties - namely elastic stiffness and fracture 
toughness, independent tuning of which would be highly 
desirable. 

While elastic responses of gels have been extensively 
studied, both in the small pQ |2] and large deformation 
regimes [H] 0, fracture studies have been up to now es- 
sentially concerned with crack nucleation and ultimate 
strength measurements [H] 0. However, trying to eluci- 
date the nature of the dissipative processes at play in frac- 
ture, which are responsible for the rate dependence of their 
strength, naturally leads to investigating the propagation 
of cracks independently from their nucleation. Tanaka et 
al [5] have performed such a study on chemical polyacry- 
lamid/ water gels. By changing the concentration of cross- 
linking agent at fixed polymer content, they found that, in 
this material, stiffness and toughness are negatively cor- 
related : as is the case for rubbers, the stiffer the gel is, 
the smaller its fracture energy. More recently, Mooney 



We report here the results of an extensive study of 
steady, strongly subsonic, mode-I (opening) crack propa- 
gation in gelatin gels. This choice was made for several 
reasons. First, due to their massive industrial use, their 
elastic properties and molecular structures have been thor- 
oughly studied. On the other hand, they can be easily cast 
into the large homogeneous samples required for fracture 
experiments. Morevover, solvent viscosity can be tuned by 
using glycerol/ water mixtures. 

We have studied the dependence of the fracture en- 
ergy Q on the crack velocity V for gels differing by their 
gelatin concentration c, glycerol content cf>, and thermal 
history, each of which is known to affect their elastic prop- 
erties. Experimental methods are described in Section 2. 
We present in Section 3. the behavior of G(V) for 3 differ- 
ent series of samples : 

A — Common c and history, variable <fi (hence solvent 
viscosity rj s ). 

B — Fixed c and </>, different histories. 

C — Common <j) and history, variable c. 

We discuss and interpret these results in Section 4. 
As already reported in JT] , the analysis of solvent effects 
(series A) leads us to propose that, in contradistinction 
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with chemical hydrogels, physical ones do not fracture by pies. So, we concluded that our preparation method does 

chain scission, but by viscous pull-out of whole gelatin not, as might have been feared, induce significant gelatin 

chains from the network via plastic yielding of the CL. hydrolysis. 

This interpretation properly accounts for the quasi-linear The pre-gel solution is poured into a mould consisting 

dependence of Q on rj s V as well as for the orders of mag- of a rectangular metal frame and two plates covered with 

nitude of its slope r = dQ/d(rj s V) and of its quasi-static Mylar films. On the longest sides of the frame, the curly 

limit Qo- We then turn toward the variations of r with part of an adhesive Velcro tape improves the gel plate 

the small strain shear modulus (i*. We find that our frac- grip. Unless otherwise specified (see Section 3.2, series B 

ture scenario, when combined with the model proposed by results), the thermal history is fixed as follows. The mould 

Joly-Duhamel et al ^21 f° r gelatin network structure and is set at 2±0.5°C for 15 h, then clamped to the mechanical 

elasticity, is compatible with the results from series B. One testing set-up and left at room temperature (19±1 °C) for 

step further, the analysis of the effect of gelatin concentra- 1 h. This waiting time ensures that variations of elastic 

tion variations (series C) leads us to invoke a concentation- moduli over the duration of the subsequent run can be 

dependent effective viscosity affecting the viscous drag on safely neglected [2] ■ The removable pieces of the mould are 

chains pulled out of the gel matrix. then taken off, leaving the 300 x 30 x 10 mm 3 gel plate fixed 

to its grips. The Mylar films are left in position to prevent 



2 Experimental methods 
2.1 Sample preparation 



solvent evaporation. They are peeled off just before the 
experiment. 

2.2 Gel characterization 



The gels are prepared by dissolving gelatin powder (type 

A from porcine skin, 300 Bloom, Sigma) in mixtures con- For each fracture experiment we prepare simultaneously 

taining a weight fraction <j) of glycerol in dcionized wa- two nominally identical samples, one of which is used 

ter, under continuous stirring for 30 min at 90°C. This to determine the elastic characteristics. For this purpose, 

temperature, higher than commonly used ones (~ 50 - with the help of the mechanical set up described below, 

60°C) has been chosen, following Ferry [T3|, so as to ob- we measure the the force-elongation response F(X) of the 

tain homogeneous pre-gel solutions even at the highest 4> plate (see Fig. QJ, up to stretching ratios A = 1.5, at the 

(60 %). A control experiment carried out with a (pure loading rate A = 1.7 10~ 2 sec -1 . 

water) /gelatin sample prepared at 60° C resulted in dif- From these data, we extract an effective small strain 

ferences of low strain moduli and r values of, respectively, shear modulus /j,* . In hydrogels, while shear stresses are 

1 % and 7 %, compatible with scatters between 90°C sam- sustained by the network, pressure is essentially borne by 
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Fig. 1. Nominal stress a = F/(eoLo) versus stretching ratio 
for a c = 10 wt%, = wt% sample plate. The dashed line is 
the extrapolation of the small strain linear response. Its slope 
is four times the effective shear modulus p* = 11 kPa (see 
text). 

the solvent. Hence, since shear moduli are typically in the 
1-10 kPa range, the gels can be considered incompress- 
ible (Poisson ratio v = 1/2), as long as no solvent draining 
occurs [T2|. So, the sound velocity relevant to define the 
subsonic regime is the transverse one c s = \J pj p, with 
p the gel mass density. For our systems, typically c s ~ 1 
m.sec -1 . Neglecting finite size effects, we assume plane 
stress uniform deformation for our plates of undeformed 
length Lq = 300 mm, width ho — 30 mm, thickness eo 
= 10 mm. In the linear regime, this assumption leads us 
to define a (necessarily somewhat overestimated) effective 
modulus as p* — j (^f) A _ ? with a = F/(eoLo) the nom- 
inal stress, A = h/ho the stretching ratio, h the stretched 
width. 

One step further, and under the conservative assump- 
tion that small strain elasticity is basically of entropic ori- 
gin, we extract a length scale characteristic of the net- 
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work as £ = (/cbT//i*) ' , which lies in the 10 nm range. 
This order of magnitude agrees with the one which can be 
evaluated from measurements of the collective diffusion 
coefficient D co u which characterizes the solvent /network 
relative motion ^3] |15| . 

For gelatin/water samples m>!> D co u ~ 10 -11 m 2 /sec, 
so that a typical time scale for draining over ~ 1cm is on 
the order of 10 7 sec, which means that macroscopic stress- 
induced draining is totally negligible here. 

As can be seen on Figure ^ beyond A values on the 
order of 1.1, the force response markedly departs from its 
small strain linear behavior. In order to calculate the me- 
chanical energy released per unit area of crack extension, 
conventionally termed energy release "rate" G, we need to 
compute the elastic energy ^-"(A) stored in the stretched 
plate. For this purpose we integrate numerically the mea- 
sured response curve. 

2.3 Fracture experiments 

The mechanical set-up is sketched on Figure |3 One of 
the grips holding the gel plate is clamped to the rigid 
external frame. The other one is attached to one end of a 
double cantilever spring of stiffness K = 43.1 x 10 3 N.m -1 . 
The other end of the spring can be displaced by a linear 
translation stage, with a 0.1/xm resolution. The deflection 
of the spring is measured by four strain gauges glued to 
the spring leaves, with a resolution of 5.10~ 2 /im. 

In most runs, the sample stiffness is much smaller than 
the spring one, and fracture occurs in the so-called fixed 
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crack linear stage 




curly tapes 



Fig. 2. Schematic representation of the mechanical setup, 
drawn around a genuine photograph of a gel plate (c = 10 
wt%, (f> — wt%), stretched to A = 1.5. Note that the crack 
propagates straight along the mid-plane. The light blue hue of 
the gel (color on line) results from Rayleigh scattering by small 
scale gel network randomness. 



Away from the sample edges, in this configuration, 
cracks run at constant velocity 1 . As expected, the free 
edges affect crack propagation up to a distance compara- 
ble with the plate width. Further data processing has been 
systematically restricted to the central region, extending 
over ~ 200 mm. In this region, we can legitimately com- 
pute the energy release rate as [I7| Q = T/(eaL ). 

Such experiments result in one run producing one sin- 
gle Q ~ V data point, hence arc very time consuming. So, 
in a second set of experiments, the stretching ratio was 
increased at the constant rate A = 1.7 10~ 2 sec -1 . This 
results in a slowly accelerating crack. We have validated 
the corresponding Q{V) data by comparison with steady 
state ones on an overlapping velocity range (see Fig. |3J) . 
The crack dynamics in this latter type of experiments can 
therefore be termed " quasi-stationary" . 



grips configuration. The stretching ratio A is computed in 
all cases by subtracting the spring deflection amplitude 
from the stage displacement. 

Before stretching, a knife cut of length 20 mm is made 
at mid-width at the upper free gel edge. In a first set 
of experiments the grips are pulled apart for 1 sec up to 
the desired amount Ah. The resulting crack advance is 
monitored by a camera with a 631 x 491 pix 2 CCD de- 
vice operating at a typical rate of 15 sec -1 . The crack tip 
position is measured with 0.5 mm resolution. The crack 
velocity V is obtained from a sliding linear regression over 
5 successive position data. 



3 Experimental results 
3.1 Solvent effects 

We summarize here the results, already reported in refer- 
ence corresponding to series A, namely gels prepared 
as described above, with gelatin concentration c — 5 wt%, 
glycerol content ranging from to 60 wt%, i.e. solvent 
viscosity rj s from 1 to 11 times that of pure water. 

1 This is true for not too small velocities, where bulk creep 
during a run is negligible. For slow cracks, with velocities below 
a few hundred /im.sec - , creep results in a measurable velocity 
drift. We only retain data out of this range. 




-xo- 
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As shown on FigureEl for all samples Q increases quasi- 
linearly with V in the explored range and, within experi- 
mental accuracy, the various curves extrapolate to a com- 
mon, ^-independent value Q(V — > 0) = Qo which yields 
an evaluated quasi-static toughness. This cannot be ac- 
cessed directly. Indeed, the above mentioned importance 
of creep in our gels leads to the well-known problems met 
when trying to define static threshold in weak solids (such 
as colloidal gels, pastes,. . . ). For this series, we find Qo ~ 
2.5 J m~ 2 , a value about 20 times smaller than a gel-air 
surface energy. 
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Fig. 4. Same data as Fig. El replotted versus r/ s V , with rj a 
the viscosity of the glycerol/water solvent. From ref. 
(reprinted from Nature Materials). 
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Indeed, the corresponding plot (Fig. 0} captures most of 
this variation. We therefore write 



g = g +rr] s v 



(i) 



The dimensionless slope r is found to be a huge number, of 
order 10 6 . In Section 4 below, we will relate the variations 
of r with those of the elastic modulus fx*. Figure |S] shows 
that, within series A, r increases with /x*. 



Fig. 3. Fracture energy release rate for gels with the same 
gelatin concentration (c = 5 wt%) and various glycerol con- 
tents (series A): (f> = wt% (circles), 20 wt% (triangles), 30 
wt% (squares), 60 wt% (diamonds). Filled symbols correspond 
to stationary cracks, open symbols to cracks accelerated in 
response to a steady increase of A. Qo = 2.5 ± 0.5 J.m~ 2 is 
the common linearly extrapolated toughness. From ref. 
(reprinted from Nature Materials). 

Moreover, the slope dQ/dV strongly increases with <fi, 
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Fig. 5. Rate sensitivity F = dQ / d(r] a V) vs. n* for the samples 
of series A. The line is the best power law fitf~ /i* 1 ' 2 . Insert 



which suggests that r] s V might be the relevant variable, shows that increasing the glycerol content stiffens the gel. 
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The quasi-scaling of Q with rj s V points toward the crit- 
ical role of polymer-solvent relative motion in the fracture 
process. In order to shed further light on this point, we 
have also performed, with the same gels, experiments in 
which a small drop of solvent is introduced into the al- 
ready moving crack opening. For such wetted cracks, as 
shown on Figure El G(V) is simply shifted downward by 
a constant amount —AQq, its slope remaining unaffected. 
The energy cost AQq ~ 2 J m~ 2 , a substantial fraction of 
Go. It clearly signals that, in the non- wetted tip case, frac- 
ture involves exposing gelatin chains to air. Such local sol- 
vent draining into the gel bulk is likely to result from the 
impossibility for our not very thin incompressible plates 
to accommodate the high strain gradients which develop 
close to the tip without being the seat of high negative 
fluid pressures. 



Fig. 6. G(V) curves for a 5 wt% gelatin gel in pure water : 
"dry" cracks opening in ambient air (upper data) and "wet" 
cracks with a drop of pure water soaking the tip. At Q too 
low for dry cracks to propagate, wet ones can still run. Linear 
fits are shown. The wet data appear merely translated towards 
lower energies. The extrapolated fracture energy for wet tips is 
Go et = 0.6±0.15 J.m~ 2 . From ref. [TTJ • (reprinted from Nature 
Materials) . 

3.2 History-controlled stiffness effects 



In a static situation, the solvent would get sucked from 
the bulk into the tip region, leading to gradual smearing 
out of the fluid pressure gradient. However, in the steadily 
moving case, the space range of this collective diffusion 
process is limited to ~ D co u/V ^S] For tip veloci- 
ties above ~ 1 mm sec" 1 , this length is smaller than the 
mesh size £, and the process is inefficient. For much slower 
cracks, it would lead to a long transient towards a lower 
apparent Qq. Trying to disentangle this from creep effects, 
which also become relevant for slow cracks, will demand a 
detailed characterization of creep which is out of the scope 
of this paper. 



The results for series A above suggest a positive correla- 
tion between the slope r and the small strain modulus fi* . 
In a second set of experiments, we have tuned \i* at two 
different gel compositions, namely cf> = 0, c = 10 and 15 
wt%. This was realized by taking advantage of the rather 
strong dependence of fi* on the temperature maintained 
during gelation, as well as on the duration of the gelation 
phase itself 52 (always chosen large enough for fi* 
variations to remain negligible during the run). This en- 
abled us to induce fi* values differing by at most a factor 
of 2. The data are shown on Figured It is seen that, for 
each c-value, again, the stiffer the gel, the tougher. Note, 
however, that r is not a function of fi* only, but also of 
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composition - a point which will be discussed in detail in 
Section 4. 
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Fig. 8. r vs. for gels from series C {4> = 0, various gelatin 
concentrations). Insert shows fi* vs. c. The full lines are the 
power law fits (see text). 



Fig. 7. r vs. [i* for gels from series B (<j> = 0, various thermal 4 DisCUSSiotl and interpretation 

histories), c = 15 wt% (full dots); c = 10 wt% (open circles). 
The curves are guide for the eye. 



4.1 A viscoplastic model of gelatin fracture 



3.3 Gelatin concentration effects 



We have investigated this last point directly by working 
with a third set of samples (series C) with the common his- 
tory described in section 2, the same solvent (pure water) 
and different values of c. As already amply documented 
HI [13 > A** increases with c (Fig(HJ) . A power law fit yields 
fi* ~ c i-64±o.2_ xhis exponent, somewhat lower than usual 
values (< 2), is close to that measured by Bot et al 6 . Fig- 
ure|S|also shows the r(fi*) data. Once more, dr/dfi* > 0. 



At first glance, as far as fracture is concerned, our gels 
share two salient features with another class of soft elastic 
materials, namely rubbers [201 |21j . In both cases : 

(1) the toughness Go is at least one order of magnitude 
larger than the energy of the surfaces created by the crack 
advance. 

(2) G increases rapidly with V in the strongly subsonic 
regime. 

Hence a first question : are the physical mechanisms 
now well established to be responsible for these two fea- 
tures in the case of rubbers also at work for our physical 
gels? 

The basic theory of rubber toughness was formulated 
by Lake and Thomas |2(J|. Fracture occurs via chain scis- 
sion : the polymer segments, of areal density S, crossing 
the fracture plane are stretched taut until they store an 
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elastic energy per monomer on the order of the covalent On the one hand, we claim that, in physical gels, frac- 

monomcr-monomer bond one, U c hain ~ a few eV. At this ture cannot process via chain scission. Indeed, the force 

stage, each of them sustains a force f chain ~ U c hain/o,, f chain defined above is more than one order of magni- 

with a a monomer size. A bond-breaking event thus cor- tude larger than that, /* ~ Ucl/q>, which can be sus- 

responds to dissipating all of the elastic energy that was tained by the H-bond stabilized cross-links. Clearly, when 

stored in the whole segment (n monomers) joining two the stored elastic energy reaches ~ Ucl per monomer, 

cross-links, ~ nV 'chain- so, gt b) ~ nU c hain^ 'j an expres- CL bonds yield, by either unzipping |27| |28| or frictional 

sion which explains the order of magnitude of Qo as well as sliding [2H]. This leads us to postulate that, in the highly 

its decrease when stiffness increases (the stiffer a rubber stressed active tip zone, the chains which cross the crack 

is, the less tough). plane creep until they are fully pulled out of the gel ma- 

The V-dependent fracture energy of rubbers is of the trix. The threshold stress at the onset of CL yielding is 

form [22] [22] (J* = f*S, with £ the areal density of crossing chains. As 

a rough estimate for this density we take £ ~ 1 /£ 2 , with 



g(^ b )(v) = g( rub) [i + $(a T v)] (2) 

where <zt is a temperature dependent WLF-like factor. V 

the above-defined estimate of the mesh size of the polymer 



This velocity dependence has been shown to result from 
bulk viscoelastic dissipation [23 [23- Due to the stress 
gradients ahead of the moving crack, which extend far be- 
yond the "active tip zone" where decohesion takes place, 
the material deforms at a strain rate which sweeps its 
whole relaxation spectrum, hence the WLF scaling fac- 
tor. 

two facts [21] : (i) linear elasticity preserves the univer- 

, _-i/9 -nii f \ ,i ii i ii Section 3.1), the plastic zone deforms under this constant 

sal r 1 stress concentration field (u) the so-called small 



network. Then, with a ~ 0.3nm, Ucl ~ O.leV, £ ~ lOnm, 

we obtain a* ~ 500kPa. 

Note that, contrary to standard conditions met with 

hard materials, here cr*//i* ~S> 1 (<~ 10 2 ), which makes the 

issue of elastic blunting raised by Hui et al [3U] certainly 

m, , n {rub) f , , . rn, u c relevant to gel fracture, 

that C/q factors out m expression (J2J) results from 

When solvent can be pumped from a wetting drop (see 



scale yielding assumption holds, namely the size of the ac- 
tive zone is negligible as compared with that of the viscous 
dissipating one. 



stress until the opening 6 C at the tip reaches the length 
of the chain - i.e. its full contour length /, since at this 
stress level it is pulled taut. This is precisely the well- 
know Dugdale model of fracture |31) . which yields, for the 
quasi-static fracture energy of wet cracks : 



We will now argue that none of these mechanisms is 
relevant in our case. Qq c1 = a* I (4) 
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From series A resuts, we estimate Q^ et w 0.6±0.15 J m~ 2 . ]TQ have shown that, for gelatin gels such as used in this 
This, together with expression enables us to get an es- work, k/£ 2 ~ 6.10 -2 . We thus estimate a V i S as resulting 
timated chain contour length I ~ 1.2/^m. With an average from the build up of the Darcy pressure over a length ~ I, 
mass M res = 80 g/mole for each of the I /a residues, this i.e. 

I S 

means a reasonable 300 kg order of magnitude estimate a vis ~ / (Vp)_D arcy ~ — — (6) 



and 



for the gelatin molar weight. 

In this picture, we interpret the shift AGo — Ga — G$ et 

as an energy cost associated with chain extraction out of S(V) ~ Go + la V i S 

I 2 

the solvent. This yields for the solvation energy per chain = Go + ot — rj s V (7) 



AGof ~ 1000 eV, i.e. ~ 10k B T per residue. 

Let us now turn to the ^-dependence of G- The tip 

wetting experiments (see Figure directly show that Go 

I 2 

and the slope r are independent : wetting shifts Go while r = a — (8) 

K 



which exhibits the observed linear variation with r) s V and 
predicts that the slope 



We found (Section 3.1) that T is of order 10 6 . With I 
as evaluated above and £ ~ 10 nm, we get from expression 
||SJ F w 2.10 5 a, which suggests that a should be of order 
1 at least. In the Dugdale model, one gets : 



leaving _T unaffected. We consider that this empirical ar- 
gument by itself rules out bulk viscoelasticity as the con- 
trolling mechanism. This appears all the more reasonable 
that rheological studies |2] show that viscous dissi- 
pation in hydrogels (loss angles typically < 0.1) is much 

S c a* 

smaller than that in rubbers. ct = — — ~ — (9) 

dact M 

For hard solids, a* is the plastic yield stress ay, al- 
ways /i. We pointed out that, for physical gels, on the 
contrary, <7*//i 3> 1. The Dugdale analysis can certainly 
not be directly used here, due to the very large defor- 
mation levels involved, hence to problems such as elas- 
tic blunting, strain-hardening and strain induced helix-coil 
transitions [35]. We were able, with the help of a hetero- 



We are therefore led to extend our fracture model to 
finite velocities. A finite V means a finite average pull-out 
velocity S — aV , where a is a geometrical factor charac- 
teristic of the shape of the Dugdale zone. Pull-out implies 
motion of the network relative to the solvent, hence a vis- 
cous contribution to the viscoplastic tip stress : 

°tip = <J* + a vis (V) (5) 



Solvent/network relative motion is diffusive |14|. which wetting experiment (pure water wetting a crack tip in a 

implies that fluid pressure gradients obey a Darcy law with glycerolled gel) reported in JT] , to obtain a direct evalu- 

an effective porosity k = rj s D co u/ /x, which can be expected ation of the size of the active zone. It yielded d ac t ~ 100 

on dimensional grounds to scale as £ 2 . Baumberger et al nm, from which we expect that a = l/d ac t ~ 10. 
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We should point out that our model for tip dissipa- As seen on Figure El the agreement with experimental 
tion (Eq. 0) is formally identical to that put forward data is quite satisfactory, bringing good support to the 
by Raphael and de Gennes [33] m the context of rubber- model. 



rubber adhesion with connector molecules. But in the gel 
case, where viscous dissipation is controlled by solvent- 
network friction, the very large compliances involved cast 
doubt on the legitimity of mathematical treatments based 
upon small opening and linear elasticity approximations 
j2U HI]. However, the possibility of accessing d act , and 
thus the fracture parameter a experimentally, together 
with the absence of substantial bulk viscoelastic dissipa- 
tion enable us to conclude that our fracture model is con- 
sistent with experiments as far as : 

- it accounts for the linear dependence of Q on r\ s V . 
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Fig. 9. Data from Fig. 0replotted versus (// 



-,2/3 



(eq. J23J) 



static toughness and the slope I '. 



Note, however, that the two data sets pertaining to 
it yields reasonable orders of magnitude for the quasi- t ^ e tw0 different gelatin concentrations do not collapse 

onto a single master curve (here a straight line). That is, 
the fracture "rate sensitivity" r does not depend on one 
single structural parameter. This remark must be consid- 
ered in the light of the finding by Joly-Duhamel et al ^21 
(hereafter abbreviated as JHAD) that, for gels of vari- 
ous gelatin concentrations, glycerol contents and thermal 
For further confirmation we now need to test the predic- histories, there is a one-to-one correspondence between 
tions of the model against the measured variations of r the storage modulus \i and the so-called helix concentra- 



4.2 Relationship between fracture and elastic 
properties 



with small strain elastic modulus /i*. 



tion Chei- This latter structural parameter, directly ob- 



T , c -j ,i i, r • t-, • i • taincd from optical activity measurements, is interpreted 

Let us first consider the results of series B, involving r j > r 

, , ,. as proportional to the length of triple-helix cross- links per 

gels with the same composition but various thermal his- t f & r i 

A j • , ibt j-i_i_ii_r i unit volume of gel. One might then be tempted to think 

tones. According to equation we predict that, for each & & r 



such set of samples, r should scale as k 1 , i.e. as 



r 



2/3 



that the modulus /i contains essentially all the mecano- 
structural information about the gel. That such is not the 
(10) case is shown by two observations : 
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(i) JHAD also show that the loss modulus fi" does not 
depend on Chei only, but also on e.g. the gelatin concenta- 
tion c. 

(ii) A non universal behavior was also found by Bot et 
al |H] for the non-linear part of the stress response in com- 
pression and in shear - a result confirmed by our own data. 

We therefore now turn to the results of series C, which 
involve gels with the same history and glycerol content 
{(f) = 0) and four different values of c. As can be seen on 
Figure HUI r/(^*) 2 / 3 definitely increases with fj,* , i.e. with 
gelatin concentration. It was shown in JHAD that, in the 
range of moduli explored here (fi > 2 kPa), gel elastic- 
ity is well described as that of a freely-hinged network of 
triple helix rods with average distance d ~ (fe^T//!) 1 / 3 , 
i.e. scaling as the mesh length scale £. This leaves the 
k^ 1 ~ ^j 2 / 3 scaling unaffected. We are thus led to at- 
tributing the residual variation of r to a concentration 
dependence of the viscosity appearing in the poroelastic 
Darcy law. We propose that this should involve, not the 
bare solvent viscosity, but an effective one 

Veff (c) = Vs e( c ) (ii) 

including possible contributions from dangling ends, 
loops attached to the network or free chains, invoked in 
JHAD and in Tanaka's study (Hj of the fracture of chemical 
gels. In view of the discussion (see Section 4.1) of the order 
of magnitude of r, clearly, 0(c) should be 0(1). 

A tentative power law fit (Figure ll(J|) yields r] e f f (c) ~ 
(/Li*)°- 75±0 03 which, combined with the /U*(c) variations 



gel as a viscoplastic disentanglement process 




5 10 15 20 25 

H* (kPa) 



Fig. 10. Data from Fig. BDreplotted as r/(//) 2/3 vs. fj,* . The 
line is the best power law fit (exponent 0.75). 

(see section 3.3), results in r] e ff(c)/r/ s ~ c . The study 
of creep viscosity in gelatin by Higgs and Ross-Murphy 3 
concluded to a c variation. However, their work was con- 
cerned with stress levels (<t//i from 2.10~ 2 to 2.1CP 1 ) con- 
siderably smaller than those relevant to the active crack 
tip zone 2 . So, though encouraging, this comparison is of 
merely indicative value. 

Finally, let us come back to the results from series A 
(same history and gelatin content, various glycerol con- 
tents (f>). A power law fit of the data shown on Figure 
yields T ~ (/i*) 1 - 2 . Here again, we must conclude that an 
increase in cj) gives rise to an increase, not only of the gel 
stiffness, but also of the effective viscosity t] e ff- Following 
2 The viscosities measured in are of order 10 s Pa sec. 
This order of magnitude, huge as compared with what we 
expect here for r/ e ff, must clearly be assigned to the stress 
range which they investigate. Indeed, far below the yield stress 
level [a <^ a*), thermally activated CL creep is necessarily 
extremely slow. 
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JHAD, an increased stiffness means an increase of Cheh References 
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which signals a change of solvent quality. In the unstressed 
gel, this most probably influences the CL average length 
as well as the helix fraction. Since changing the Flory in- 
teraction parameter shifts helix-coil transitions, it is likely 
to also affect the structural changes shown by Courty et 
al |32| to result in large variations of optical activity in 
the large strain regime. We expect the value of rj e f / to be 
sensitive to these structural modifications. 

In conclusion, we contend here that fracture of chem- 
ical and physical gels is controlled by different mecha- 
nisms : 

- stretched chain scission (chemical gels). 

- viscoplastic cross-link yield leading to chain pull-out 
(physical gels). 

Of course, the model formulated here should be tested 
more completely by studying crack tip dynamics in other 
physical hydrogels involving CL with different structures, 
such as ionically cross-linked alginates. More work will 
also be needed along two directions : (a) characterization 
of creep dynamics at larger stress levels than those used in 
reference [3], and of its dependence on solvent viscosity; 
(b) more detailed study of slow crack motion, aimed at 
improving the reliability of (Jo-determinations as well as 
at testing possible effects of bulk poroelasticity. 

We are gratelul to C.Y. Hui for an enlightening discussion. We 
thank L. Legrand for his contribution to the analysis of the gel 
light-scattering properties. 
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